Active Model-Based Hysteresis Compensation and Tracking Control of Pneumatic Artificial Muscle

The hysteretic nonlinearity of pneumatic artificial muscle (PAM) is the main factor that degrades its tracking accuracy. This paper proposes an efficient hysteresis compensation method based on the active modeling control (AMC). Firstly, the Bouc–Wen model is adopted as the reference model to describe the hysteresis of the PAM. Secondly, the modeling errors are introduced into the reference model, and the unscented Kalman filter is used to estimate the state of the system and the modeling errors. Finally, a hysteresis compensation strategy is designed based on AMC. The compensation performances of the nominal controller with without AMC were experimentally tested on a PAM. The experimental results show that the proposed controller is more robust when tracking different types of trajectories. In the transient, both the overshoot and oscillation can be successfully attenuated, and fast convergence is achieved. In the steady-state, the proposed controller is more robust against external disturbances and measurement noise. The proposed controller is effective and robust in hysteresis compensation, thus improving the tracking performance of the PAM.


Introduction
Pneumatic artificial muscle (PAM) is a widely-utilized bionic flexible actuator. PAM is mainly composed of a hollow rubber tube, a fiber woven mesh and metal connectors. PAM-driven robots are more flexible than hydraulic or motor-driven robots. In addition, PAM features a higher power-to-mass ratio and good compliance with the human body. Therefore, it has been widely used in the medical and rehabilitation fields. PAM is used in the exoskeleton system, and a terminal sliding mode control was adopted to higher trajectory tracking accuracy [1]. Although PAMs feature good application prospects, the inherent nonlinear hysteresis characteristic has become one of the main obstacles affecting the precise trajectory tracking control of PAM-driven robots. Currently, hysteresis modeling and compensation have gradually become significant areas of research focus [2][3][4][5][6].
From the mathematical point of view, dynamic models for PAMs can be briefly divided into physics-based and phenomenon-based models [7]. The physics-based model mainly analyzes the geometric relationship and material properties of the PAM. For example, C. Kothera et al. considered the thickness of the rubber tube and used the force balance method to derive a static model of a PAM [8]. Because the physics-based model involves very complex mathematical derivation and a large number of model parameters, it creates great difficulties for the design of the controller and adds a very large calculation cost.
Alternately, models based on phenomenology are obtained through experimental data, making the description of the model simpler. For example, in the research work of R. Colbrunn et al., the PAM is regarded as a structure composed of a damping unit, a spring unit and a Coulomb friction unit in parallel [9]. In previous studies, many When the inner pressure decreases, the PAM gradually returns to its initial state under the action of the outer braided mesh. Due to the friction between the rubber tube and the braided mesh, there is strong hysteretic nonlinearity between the inner pressure and the output displacement and force.
For instance, for the PAM utilized in Figure 1a, the control voltage to the proportional pressure regulator is adopted as the input and the displacement of the PAM is adopted as the output. Figure 1b shows the measured input-output hysteresis loops of the PAM. It can be observed that the hysteresis loop of the PAM is not symmetric, and the loop shape is influenced by the frequency of the control input, i.e., the rate-dependence. These characteristics increase the difficulty of precisely modeling the hysteresis.

Bouc-Wen Hysteresis Model
In this paper, the Bouc-Wen model is used to describe the hysteresis characteristics of the PAM. The Bouc-Wen model was originally proposed for nonlinear vibration mechanics [25], while currently it is a popular hysteresis model. In the Bouc-Wen model, a state variable h is used to characterize the hysteresis relationship of the system. The relationship between the input and the state variable h can be described in the following form: where u is a generalized input, h is the hysteresis state variable, α, β and γ are the gains controlling the shape of the hysteresis loop and m controls the smoothness of the transition from elastic to plastic response. In applications, m can be set to 1 so as to simplify the model structure. This is also adopted in this paper. For the PAM investigated in this paper, the control voltage applied to the proportional pressure regulator valve is the input of the overall system. As a result, the control voltage is adopted as the input in the Bouc-Wen model. Therefore, the Bouc-Wen model for the PAM can be expressed as follows: For instance, for the PAM utilized in Figure 1a, the control voltage to the proportional pressure regulator is adopted as the input and the displacement of the PAM is adopted as the output. Figure 1b shows the measured input-output hysteresis loops of the PAM. It can be observed that the hysteresis loop of the PAM is not symmetric, and the loop shape is influenced by the frequency of the control input, i.e., the rate-dependence. These characteristics increase the difficulty of precisely modeling the hysteresis.

Bouc-Wen Hysteresis Model
In this paper, the Bouc-Wen model is used to describe the hysteresis characteristics of the PAM. The Bouc-Wen model was originally proposed for nonlinear vibration mechanics [25], while currently it is a popular hysteresis model. In the Bouc-Wen model, a state variable h is used to characterize the hysteresis relationship of the system. The relationship between the input and the state variable h can be described in the following form: where u is a generalized input, h is the hysteresis state variable, α, β and γ are the gains controlling the shape of the hysteresis loop and m controls the smoothness of the transition from elastic to plastic response. In applications, m can be set to 1 so as to simplify the model structure. This is also adopted in this paper. For the PAM investigated in this paper, the control voltage applied to the proportional pressure regulator valve is the input of the overall system. As a result, the control voltage is adopted as the input in the Bouc-Wen model. Therefore, the Bouc-Wen model for the PAM can be expressed as follows: where Y(t) is the PAM's output displacement, u is the control voltage, ρ ≤ 1 is the weight coefficient, and k is the stiffness coefficient of the PAM [26]. The schematic diagram of the hysteresis model is shown in Figure 2, where the parameters of the Bouc-Wen model are shaded in grey background. where Y(t) is the PAM's output displacement, u is the control voltage, ρ ≤ 1 is the weight coefficient, and k is the stiffness coefficient of the PAM [26]. The schematic diagram of the hysteresis model is shown in Figure 2, where the parameters of the Bouc-Wen model are shaded in grey background.

Active Modeling for PAM
In this paper, the Bouc-Wen model defined in (2) is adopted as the reference model in describing the PEA's hysteresis. Compared with the operator-based hysteresis models, the Bouc-Wen model requires fewer parameters, which is beneficial to the controller design. However, the modeling accuracy of the Bouc-Wen model is limited, especially when it is used to describe the hysteresis characteristics at different pressure ranges or at different frequencies. In order to improve the modeling accuracy, we introduce modeling errors into the reference model. The state and modeling error of the reference model are taken as the extended state vector and a UKF is used to estimate it in real time. Equation (2) can be rewritten as: where U(t) and V(t) are the process noise and the observation noise of the system. Due to interference and unmodeled uncertainties, there are modeling errors between the actual system and the Bouc-Wen model. We use fh(t) to depict this modeling error, and thus the model of the actual system can be described as: where  ( ) h t and  ( ) Y t are the state of the actual dynamics of the PAM and the output displacement of the system. Currently, closed-loop systems often run at a high sampling rate, typically 1 kHz or above. For this high sampling rate, the model error fh(t) can be treated as a slowly changing parameter and it can be approximated as: where h(t) is the state of the reference model in (2) and p(t) is assumed to be the process noise actuating the model errors. Subsequently, the discrete description of (4) is obtained:

Active Modeling for PAM
In this paper, the Bouc-Wen model defined in (2) is adopted as the reference model in describing the PEA's hysteresis. Compared with the operator-based hysteresis models, the Bouc-Wen model requires fewer parameters, which is beneficial to the controller design. However, the modeling accuracy of the Bouc-Wen model is limited, especially when it is used to describe the hysteresis characteristics at different pressure ranges or at different frequencies. In order to improve the modeling accuracy, we introduce modeling errors into the reference model. The state and modeling error of the reference model are taken as the extended state vector and a UKF is used to estimate it in real time. Equation (2) can be rewritten as: where U(t) and V(t) are the process noise and the observation noise of the system. Due to interference and unmodeled uncertainties, there are modeling errors between the actual system and the Bouc-Wen model. We use f h (t) to depict this modeling error, and thus the model of the actual system can be described as: where h(t) and Y(t) are the state of the actual dynamics of the PAM and the output displacement of the system.
Currently, closed-loop systems often run at a high sampling rate, typically 1 kHz or above. For this high sampling rate, the model error f h (t) can be treated as a slowly changing parameter and it can be approximated as: where h(t) is the state of the reference model in (2) and p(t) is assumed to be the process noise actuating the model errors. Subsequently, the discrete description of (4) is obtained: where u δ(k) = (u (k) − u (k−1) )/T s is the discrete expression of the derivative of the input voltage, T s is the sample time, h k , h k , u k , f h(k) and Y k are the discrete expressions of h(t), h(t), u(t), f h (t) and Y(t), respectively and U k and V k are the sampling value of U(t) and V(t), respectively. In this paper, the "joint estimation" technique is adopted [27]. The discrete extended state is defined as follows: In order to use UKF, we rewrite (7) as: where F and H are the state update function and the measurement function described as: Due to the nonlinear characteristics of the Bouc-Wen model in (2), UKF can be used as the estimator. UKF is a well-developed filter widely used for estimation, as in our previous work [22]. Combined with the Bouc-Wen model, UKF includes the following steps.
The first step is to calculate the sampling points based on the estimated mean and covariance: where n is the state dimension, λ is a constant and χ k−1 is the matrix after sampling expansion, calculating from the estimated mean and covariance calculated on the state vector X k−1 . The second step is the prediction: where χ k|k−1 is obtained by the nonlinear conversion of sampling points χ k−1 , Q k is the covariance matrix of process noise and the process noise is added to the covariance prediction step to obtain the weighted covariance matrix in the state space P X(k|k−1) . The third step is the update: where R k is the covariance matrix of the measurement noise, Υ χ(k|k−1) is the measurement prediction of χ k−1 from time k − 1 to time k, K k is the Kalman gain, X k|k is the estimation of state vector based on confidence field χ k|k−1 and P k|k is the update of confidence matrix P k−1|k−1 . W c (i) and W m (i) are described below: where λ = σ 2 (n + κ), κ is the scale factor and its value only needs to ensure that the covariance matrix is non-negatively definite, σ controls the range of sampling points distribution and the tuning of ϑ can improve the approximate accuracy of the covariance. The tuning of the parameters is straightforward. The range of σ is 10 −4 ≤ σ ≤ 1. A larger σ value helps to increase the noise suppression, whereas the robustness might become weaker. As a result, σ is usually set to be a smaller value. For Gaussian distribution, ϑ can usually be set to 2.

Active-Model-Based Control Strategy
The closed-loop controller proposed in this paper can be divided into two parts. As shown in Figure 3, the first part is the nominal controller. In this paper, the widely-utilized PID controller is selected as the nominal controller. The other part is the state estimation of the system through UKF. For the proposed AMC-based controller, the estimated state of the system is used to compensate the error of the nominal controller via the compensation strategy unit. In this way, the uncertainties of the nominal controller are compensated, and the control performance of the overall system can be improved. In this control strategy, the requirement on the accuracy of the reference model does not need to be too high, because the UKF and compensation strategy will estimate and compensate for the model uncertainty in real time.  The block diagram of the control strategy is shown in Figure 3. The total control voltage of the controller is where u is the final control voltage applied to the PAM, un and uc are the control voltages generated by the nominal controller and the compensation strategy, respectively. The control objective is expressed as: where  k Y and Yd(k) are the actual and desired trajectories, respectively, and the estimated tracking error is defined as In this paper, the PI controller is adopted as the nominal controller, which can be expressed in the following form: The block diagram of the control strategy is shown in Figure 3. The total control voltage of the controller is where u is the final control voltage applied to the PAM, u n and u c are the control voltages generated by the nominal controller and the compensation strategy, respectively. The control objective is expressed as: where Y k and Y d(k) are the actual and desired trajectories, respectively, and the estimated tracking error is defined as e k = Y d(k) − Y k . In this paper, the PI controller is adopted as the nominal controller, which can be expressed in the following form: where k p and k i are the proportional and integral gains. Based on the control voltage of the nominal controller and the result of the UKF estimation, the estimated value ofˆ Y k+1 can be obtained, which can be expressed as follows: In this paper, after getting the prediction error, we define the compensation strategy as follows where k c a proportional gain. Thus, the final control input u is obtained:

Experimental Setup
In order to verify the actual performance of the proposed control method, multiple sets of trajectory tracking experiments are implemented on an in-house built testbench for PAMs. The experimental setup is shown in Figure 4

Experimental Setup
In order to verify the actual performance of the proposed control method, multiple sets of trajectory tracking experiments are implemented on an in-house built testbench for PAMs. The experimental setup is shown in Figure 4, where different PAMs are installed on a linear rail. The Novotechnik displacement sensor is used to measure the displacement of the PAM with an accuracy of 0.001 mm. In this paper, the nominal stroke of the selected PAM (model DMSP-20-200N from Festo) is 40 mm at a maximum pressure of 0.8 MPa. The proportional pressure regulator valve (model VPPM-6L-L from Festo) is used to control the inner pressure of the PAM. The data acquisition and closed-loop control are implemented on a real-time target machine (Model Mobile from Speedgoat) at a sampling rate of 1 kHz.

Identification of the Bouc-Wen Model
As shown in (2), there are six parameters to be identified in the Bouc-Wen model, i.e., k, ρ, b, α, β and γ. Genetic algorithm (GA) is adopted as the identification algorithm. In the open-loop condition, a sinusoidal signal of u0(t) = 2.9 sin (0.2πt) + 0.2 is applied to the pressure regulator valve as the control voltage. The displacement of the PAM is measured and used for parameter identification. We use the mean square error (MSE) in the

Identification of the Bouc-Wen Model
As shown in (2), there are six parameters to be identified in the Bouc-Wen model, i.e., k, ρ, b, α, β and γ. Genetic algorithm (GA) is adopted as the identification algorithm. In the open-loop condition, a sinusoidal signal of u 0 (t) = 2.9 sin (0.2πt) + 0.2 is applied to the pressure regulator valve as the control voltage. The displacement of the PAM is measured and used for parameter identification. We use the mean square error (MSE) in the following form as the fitness function for the GA: where N is the number of sampling point, i is the sampling index, Y actual is the measured displacement and Y model is the Bouc-Wen model output.
The identified results are k = 6.835, ρ = 0.987, b = −0.569, α = −0.188, β = 1.106 and γ = 0.838. The measured displacement and the model output are plotted in Figure 5a. The model output can match the measurement, whereas the modeling accuracy is not very high.
As listed in the first row of Table 1, the RMSE values are calculated to be 2.012 mm, 3.687 mm and 2.282 mm for u0, u1 and u2, respectively.

Active Model Error Estimation
We define the expansion vector in (8) using the actual state of the system and the modeling error. Next, the expansion vector is introduced into the reference model. We use the UKF defined by (9) to estimate the expansion vector in real time, that is, to estimate the actual state and the modeling error of the system at the same time. The modeling accuracy of the active modeling is also investigated and plotted in Figure 5 as a comparison. Compared with the reference model, the system output estimated by the UKF almost coincides with the measured displacement of the system. Similarly, the RMSE is also adopted to evaluate the modeling accuracy of the active model. The RMSE values for u0(t), u1(t), and u2(t) are calculated to be 1.4822 × 10 −4 mm, 1.379 × 10 −4 mm and 1.3974 × 10 −4 mm, respectively. Because the resolution of the displacement sensor is 0.01 mm, such small magnitude errors can be approximated to zero. The RMSE of the active model is also listed in the second row of Table 1.  Based on the above experimental results, we can conclude that the reference model, i.e., the Bouc-Wen model, can better account for the frequency variation than the amplitude variation. However, it is difficult to maintain high modeling accuracy. By contrast, with the help of the expansion vector, UKF is effective in estimating the modeling errors. The modeling accuracy can be significantly improved using the proposed active model.

Extended State Observer Based Controller for Comparison
The AMC-based PI controller was established in Section IV. In order to test the performance of the proposed controller, trajectory tracking performances with and without active modeling were experimentally compared. For the purpose of comparison, an extended state observer (ESO) was also integrated with the nominal PID controller, i.e., PID + ESO. In order to use the state observer, the state equation of the system is expressed as follows: In order to verify the modeling accuracy of the identified Bouc-Wen model, a further two sinusoidal control signals with different amplitudes and frequencies are used to excite the PAM, i.e., u 1 (t) = 1.4 sin (0.2πt) + 0.2 and u 2 (t) = 2.9 sin (0.1πt) + 0.2. The measured displacements and model outputs are presented in Figure 5b,c. In Figure 5b, the amplitude of u 1 is different from u 0 , leading to a significant drop on the modeling accuracy of the Bouc-Wen model. In Figure 5c, the frequency of u 2 is different from u 0 . In this case, the modeling accuracy is comparable to u 0 . The following root mean square error (RMSE) is used to quantitatively evaluate the modeling errors: As listed in the first row of Table 1, the RMSE values are calculated to be 2.012 mm, 3.687 mm and 2.282 mm for u 0 , u 1 and u 2 , respectively.

Active Model Error Estimation
We define the expansion vector in (8) using the actual state of the system and the modeling error. Next, the expansion vector is introduced into the reference model. We use the UKF defined by (9) to estimate the expansion vector in real time, that is, to estimate the actual state and the modeling error of the system at the same time. The modeling accuracy of the active modeling is also investigated and plotted in Figure 5 as a comparison. Compared with the reference model, the system output estimated by the UKF almost coincides with the measured displacement of the system. Similarly, the RMSE is also adopted to evaluate the modeling accuracy of the active model. The RMSE values for u 0 (t), u 1 (t), and u 2 (t) are calculated to be 1.4822 × 10 −4 mm, 1.379 × 10 −4 mm and 1.3974 × 10 −4 mm, respectively. Because the resolution of the displacement sensor is 0.01 mm, such small magnitude errors can be approximated to zero. The RMSE of the active model is also listed in the second row of Table 1.
Based on the above experimental results, we can conclude that the reference model, i.e., the Bouc-Wen model, can better account for the frequency variation than the amplitude variation. However, it is difficult to maintain high modeling accuracy. By contrast, with the help of the expansion vector, UKF is effective in estimating the modeling errors. The modeling accuracy can be significantly improved using the proposed active model.

Extended State Observer Based Controller for Comparison
The AMC-based PI controller was established in Section IV. In order to test the performance of the proposed controller, trajectory tracking performances with and without active modeling were experimentally compared. For the purpose of comparison, an extended state observer (ESO) was also integrated with the nominal PID controller, i.e., PID + ESO. In order to use the state observer, the state equation of the system is expressed as follows: where u ESO is the control input, y is the actual measured value of the system, x 1 and x 2 are the state variables of the system and τ is treated as a total distribution including the unmodeled nonlinearities and external disturbance. The value τ can be extended as an additional state variable, i.e., x 3 = τ. The derivative of x 3 is represented by p 0 . Equation (22) can be rewritten as: The extended system model described in (23) is observable. A ESO proposed in [29] is used to estimate system states, which can be constructed as follows:  (24) is combined with the PID controller, the final control voltage is expressed as: where u n is the control voltage of the nominal PID controller andx 3 is the estimated state of the total disturbance of the system. The block diagram of PID + ESO is shown in Figure 6.   (24) is combined with the PID controller, the final control voltage is expressed as: where un is the control voltage of the nominal PID controller and 3 x is the estimated state of the total disturbance of the system. The block diagram of PID + ESO is shown in Figure  6.

Experimental Results
The triangular trajectory is a widely utilized continuous signal. In this section, a 0.1 Hz, 0.2-30 mm triangular trajectory is adopted as the desired trajectory so as to tune the parameters of the controllers. The gains of the nominal PI controller are tuned to be kp = 0.1 and ki = 0.15 for a better balance between the transient and steady-state tracking performances. The following parameters are initialized for the UKF in the active model:

Experimental Results
The triangular trajectory is a widely utilized continuous signal. In this section, a 0.1 Hz, 0.2-30 mm triangular trajectory is adopted as the desired trajectory so as to tune the parameters of the controllers. The gains of the nominal PI controller are tuned to be k p = 0.1 and k i = 0.15 for a better balance between the transient and steady-state tracking performances. The following parameters are initialized for the UKF in the active model: The gains of the compensation strategy are set to k p = 0.1 and k i = 0.15, k c = 0.06. For the PID + ESO, the same nominal PI controller is adopted, and the parameters in ESO are tuned to β 1 = β 2 = β 3 = 100, δ = 500.
The tracking performance of the nominal PID controller, PID + AMC and PID + ESO are shown in Figure 7a. For the open-loop system, the tracking performance is poor due to the influence of the hysteresis nonlinearity. For the nominal PI controller, the PAM can follow the desired trajectory well, except for the slight overshoots at the turning points of the trajectory. For the PID + AMC and PID + ESO controllers, the overshoots at the turning points can be suppressed. Because this trajectory is a slow-varying trajectory, the steady-state performances of the controllers are comparable. Similar to the evaluation of the modeling accuracy, the RMSE is utilized to quantitively evaluate the tracking performance of the controllers. The RMSEs of the PID, PID + ESO, and PID + AMC controllers are calculated to be 0.3461 mm, 0.3243 mm and 0.3209 mm, respectively.
The hysteresis compensation performances are illustrated in Figure 7b. Compared with the open-loop system, the hysteresis loop of the closed-loop system stays close to the 45 • line. This demonstrates that all the controllers can successfully compensate the PAM's hysteresis. Figure 7c shows the control voltages from the nominal controller and the compensation unit. It can be found that the with the help of active modeling and compensation, the control input at the turning points can be compensated.
Based on the above results, we can conclude that the parameters of all the controller have been well tuned. In order to test the robustness of the controllers, all the parameters are fixed, and tracking experiments of different trajectories are performed. Sinusoidal trajectory with varying amplitude is a very good choice for evaluating the hysteresis compensation performance as it can show both the major and minor loops of the PAM's hysteresis.
The tracking results of a 0.05 Hz sinusoidal trajectory with descending amplitude are shown in Figure 8a. Similar to the performance of the triangular trajectory, the tracking performance of the closed-loop system can be significantly improved by the controllers. the trajectory. For the PID + AMC and PID + ESO controllers, the overshoots at the turning points can be suppressed. Because this trajectory is a slow-varying trajectory, the steadystate performances of the controllers are comparable. Similar to the evaluation of the modeling accuracy, the RMSE is utilized to quantitively evaluate the tracking performance of the controllers. The RMSEs of the PID, PID + ESO, and PID + AMC controllers are calculated to be 0.3461 mm, 0.3243 mm and 0.3209 mm, respectively.
The hysteresis compensation performances are illustrated in Figure 7b. Compared with the open-loop system, the hysteresis loop of the closed-loop system stays close to the 45° line. This demonstrates that all the controllers can successfully compensate the PAM's hysteresis. Figure 7c shows the control voltages from the nominal controller and the compensation unit. It can be found that the with the help of active modeling and compensation, the control input at the turning points can be compensated. Based on the above results, we can conclude that the parameters of all the controller have been well tuned. In order to test the robustness of the controllers, all the parameters are fixed, and tracking experiments of different trajectories are performed. Sinusoidal trajectory with varying amplitude is a very good choice for evaluating the hysteresis compensation performance as it can show both the major and minor loops of the PAM's hysteresis.
The tracking results of a 0.05 Hz sinusoidal trajectory with descending amplitude are shown in Figure 8a. Similar to the performance of the triangular trajectory, the tracking performance of the closed-loop system can be significantly improved by the controllers. However, for the nominal controller, there are obvious overshoot and oscillations in the  For the hysteresis compensation, as shown in Figure 8b, the nominal controller is efficient at hysteresis compensation, whereas the obvious overshoot and oscillations are not desired in real implementations. It can also be clearly observed in Figure 8b that the proposed PID + AMC controller can improve the transient performance. Figure 8c shows the effect of the active modeling. In the transient, based on the estimated tracking error, the compensation unit can immediately compensate for the nominal controller. Experiments were further conducted on tracking another 0.1 Hz sinusoidal trajectory with descending amplitude to test the performances of the controllers for trajectories at different frequencies. The experimental results are provided in Figure 9. Similar to the results in Figure 8, the nominal controller can compensate the hysteresis and the transient tracking performances can be further improved with active modeling.
The statistics on the tracking errors of the controllers in the above experiments are summarized in Table 2. In each case, the proposed controller can further improve the performance of the nominal controller.

Open-Loop
PID PID + ESO PID + AMC For the hysteresis compensation, as shown in Figure 8b, the nominal controller is efficient at hysteresis compensation, whereas the obvious overshoot and oscillations are not desired in real implementations. It can also be clearly observed in Figure 8b that the proposed PID + AMC controller can improve the transient performance. Figure 8c shows the effect of the active modeling. In the transient, based on the estimated tracking error, the compensation unit can immediately compensate for the nominal controller.
Experiments were further conducted on tracking another 0.1 Hz sinusoidal trajectory with descending amplitude to test the performances of the controllers for trajectories at different frequencies. The experimental results are provided in Figure 9. Similar to the results in Figure 8, the nominal controller can compensate the hysteresis and the transient tracking performances can be further improved with active modeling. It must be pointed out that triangular and sinusoidal trajectories are continuous, i.e., no sudden change occurs in the trajectory. In this case, the hysteresis compensation is relatively simple. In order to test the performance of the controllers in tracking discontinuous trajectories, square wave and sawtooth trajectories are adopted as the desired trajectories. Generally, the controller parameters need to be further tuned for a better performance. In order to test the robustness of the proposed controller, the same controller parameters are inherited. The experimental results are shown in Figure 10. Figure 10a shows the tracking performance of the square wave trajectory. For the nominal controller, the transient performance is very poor. For instance, on the forward motion of the PAM, the nominal controller is fast, whereas large overshoot (~17.3%) is observed before the system enters into the steady-state. The settling time was found to be approximately 0.8 s. PID + ESO reduces the overshoot to approximately 5.7%. The proposed PID + AMC can track the square wave signal with negligible overshoot and it can converge within 0.5 s. On the backward motion of the PAM, the performances of all the controllers are comparable, whereas the proposed PID + AMC can still converge quickly without overshoot. The sawtooth trajectory is another popular discontinuous trajectory. As shown in Figure 10b, the nominal controller responds slowly at the turning point of the trajectory. Similar to Figure 10a, the proposed PID + AMC controller offers the best transient performance. For discontinuous trajectories, the sudden change in the trajectory might degrade the transient performance of the nominal controller. It can be seen from the compensation voltage in Figure 10 that the compensation effect of the proposed controller is more obvious at the sudden change of the trajectory. The statistics on the tracking errors of the controllers in the above experiments are summarized in Table 2. In each case, the proposed controller can further improve the performance of the nominal controller. It must be pointed out that triangular and sinusoidal trajectories are continuous, i.e., no sudden change occurs in the trajectory. In this case, the hysteresis compensation is relatively simple. In order to test the performance of the controllers in tracking discontinuous trajectories, square wave and sawtooth trajectories are adopted as the desired trajectories. Generally, the controller parameters need to be further tuned for a better performance. In order to test the robustness of the proposed controller, the same controller parameters are inherited. The experimental results are shown in Figure 10. Figure 10a shows the tracking performance of the square wave trajectory. For the nominal controller, the transient performance is very poor. For instance, on the forward motion of the PAM, the nominal controller is fast, whereas large overshoot (~17.3%) is observed before the system enters into the steady-state. The settling time was found to be approximately 0.8 s. PID + ESO reduces the overshoot to approximately 5.7%. The proposed PID + AMC can track the square wave signal with negligible overshoot and it can converge within 0.5 s. On the backward motion of the PAM, the performances of all the controllers are comparable, whereas the proposed PID + AMC can still converge quickly without overshoot. The sawtooth trajectory is another popular discontinuous trajectory. As shown in Figure 10b, the nominal controller responds slowly at the turning point of the trajectory. Similar to Figure 10a, the proposed PID + AMC controller offers the best transient performance. For discontinuous trajectories, the sudden change in the trajectory might degrade the transient performance of the nominal controller. It can be seen from the compensation voltage in Figure 10 that the compensation effect of the proposed controller is more obvious at the sudden change of the trajectory. Further, a payload of approximately 40 N was used to test the robustness of the controllers against external disturbances. It is exerted to the hook when the closed-loop system enters into the steady-state. In this case, it can be regarded as a constant external disturbance. The experimental results are shown in Figure 11. It can be observed that the proposed PID + AMC can quickly converge to its steady-state value, showing excellent disturbance rejection capability. By contrast, it takes longer for PID and PID + ESO to converge. As a result, the proposed controller is more robust against external disturbances. In the above experimental results, the most important improvement to the proposed controller appears in the transient performance. For the nominal controller, it is difficult to suppress the overshoot and oscillations in tracking different trajectories while maintaining a fast response capability. For the proposed AMC-based controller, due to the high-precision estimation of the modeling error, the tracking performance of the nominal Further, a payload of approximately 40 N was used to test the robustness of the controllers against external disturbances. It is exerted to the hook when the closed-loop system enters into the steady-state. In this case, it can be regarded as a constant external disturbance. The experimental results are shown in Figure 11. It can be observed that the proposed PID + AMC can quickly converge to its steady-state value, showing excellent disturbance rejection capability. By contrast, it takes longer for PID and PID + ESO to converge. As a result, the proposed controller is more robust against external disturbances. Further, a payload of approximately 40 N was used to test the robustness of the controllers against external disturbances. It is exerted to the hook when the closed-loop system enters into the steady-state. In this case, it can be regarded as a constant external disturbance. The experimental results are shown in Figure 11. It can be observed that the proposed PID + AMC can quickly converge to its steady-state value, showing excellent disturbance rejection capability. By contrast, it takes longer for PID and PID + ESO to converge. As a result, the proposed controller is more robust against external disturbances. In the above experimental results, the most important improvement to the proposed controller appears in the transient performance. For the nominal controller, it is difficult to suppress the overshoot and oscillations in tracking different trajectories while maintaining a fast response capability. For the proposed AMC-based controller, due to the high-precision estimation of the modeling error, the tracking performance of the nominal In the above experimental results, the most important improvement to the proposed controller appears in the transient performance. For the nominal controller, it is difficult to suppress the overshoot and oscillations in tracking different trajectories while maintaining a fast response capability. For the proposed AMC-based controller, due to the high-precision estimation of the modeling error, the tracking performance of the nominal controller can be further improved. During the transient, no obvious overshoot and oscillation exists, and the controller converges faster. We can conclude that high hysteresis compensation efficiency and high robustness are achieved using the proposed controller.

Conclusions
How to model and compensate the strong hysteresis nonlinearities of the PAM has become one of the main obstacles to the high-precision motion control of PAMs. This paper presents an efficient hysteresis compensation strategy, in which the Bouc-Wen model is adopted as the reference model and a UKF is utilized to estimate the state and the modeling error of the reference model. The estimated error is then utilized as a compensation to the nominal controller so as to improve the performance of the nominal controller.
Currently, the sampling rate of the closed-loop system can be set high, e.g., 1 kHz or above. In this case, in the steady-state, the estimated tracking error at the next time interval is close to the current tracking error, especially when the PAM is tracking slow and continuous trajectory. However, if the PAM is confronted with a discontinuous trajectory or disturbance, the discrepancy between the estimated tracking error and the current measurement is large. This might lead to poor transient performance. The experimental results demonstrate that for continuous and discontinuous trajectories, the transient performances can be improved by the proposed controller. The proposed controller can successfully eliminate the overshoot and oscillation of the nominal controller and fast convergence is achieved. For the steady-state performance, the proposed controller can better suppress the influence of the measurement noise. Therefore, the effectiveness of the proposed AMC-based hysteresis compensation strategy is verified.